Quantum Hall fractions in ultracold fermionic vapors 
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We study the quantum Hall states that appear in the dilute limit of rotating ultracold fermionic 
gases when a single hyperfine species is present. We show that the p-wave scattering translates into 
a pure hard-core interaction in the lowest Landau level. The Laughlin wavefunction is then the exact 
ground state at filling fraction u=l/S. We give estimates of some of the gaps of the incompressible 
liquids for v = p/(2p ± 1). We estimate the mass of the composite fermions at v — 1/2. The width 
of the quantum Hall plateaus is discussed by considering the equation of state of the system. 



PACS numbers: 03.75.Ss, 05.30.Fk, 73.43.-f, 73.43.Cd 

The preparation and manipulation of ultracold atomic 
gases has led to many interesting developments in the 
study of quantum fluids undergoing fast rotation 0, Q . 
The Bose-Einstein condensates can be set in various rota- 
tion regimes with the characteristic response of a super- 
fluid. The condensate does not acquire angular momen- 
tum below some velocity threshold. Then there is nu- 
cleation of one vortex and with increasing velocity more 
and more vortices are created. They have been observed 
forming the Abrikosov triangular lattice 0, 0. When 
the rotation frequency reaches the trapping frequency 
in the radial plane, it has been predicted that quantum 
Hall fractional states should become ground states of the 
system [a, la 11 la if the gas enters a two-dimensional 
regime. Trapped Fermi gases may also exhibit superflu- 
idity if they undergo BCS pairing condensation. If the 
pairing strength is varied, it may be possible to observe 
the crossover from molecular condensation at strong cou- 
pling to BCS phase transition at weak coupling. A possi- 
ble signature of the superfluid paired phase of fermions is 
the peculiar response to stirring, leading again to vortex 
formation. In the fast rotation limit, it is thus natural 
to ask if there is formation of fractional quantum Hall 
states as in the Bose case and what are their properties. 

In this Letter we investigate the fractional quantum 
Hall effect (FQHE) appearing in atomic vapor made of 
a single hyperfine species of fermions. We show that the 
p-wave scattering between fermions can lead to the for- 
mation of the Jain principal sequence of FQHE fractions 
v = p/(2p± 1), in addition to the celebrated Laughlin 
wavefunction at v = 1/3, as well as a Fermi sea of com- 
posite fermions for half-filling of the lowest Landau level 
(LLL) . We give estimates of the gaps for the incompress- 
ible fluids governed by the p-wave scattering length and 
of the mass of the composite fermions. The equation of 
state of the system seen as the angular momentum of the 
ground state as a function of the rotation frequency dis- 
plays plateaus corresponding to the FQHE fluids. Their 
widths can be estimated by taking into account the nu- 
cleation of quasiparticles. 

We consider a gas of fermionic atoms and suppose that 
they are set in rotation for example by a stirring external 



potential @ that can be applied for some time to transfer 
angular momentum to the gas and then is removed. We 
are then left with a rotating cloud and we assume that it 
attains thermal equilibrium in the rotating frame. If TL 
stands for the Hamiltonian in the laboratory frame then 
it becomes TLr =T~t — loL z in the rotating frame where lo 
is the rotation frequency and L z the angular momentum 
along the rotation axis. The Hamiltonian describing N 
particles of mass m in this frame can be written as : 
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where the xy trap frequency is u>q, the axial frequency 
is lo z and the angular velocity vector is ujz. For u> close 
to wo, the physics is that of charge-e particles in a mag- 
netic field B = (2mw/e)z, corresponding to a magnetic 
length £ = (2mw). We assume the existence of a 
two-dimensional (2D) regime in which the wavefunction 
along the z-axis is the ground state of the z-axis harmonic 
potential. 

If we consider a single hyperfine species of fermions, 
then the s-wave scattering is forbidden by the Pauli prin- 
ciple. The next allowed partial wave, the p-wave, leads to 
much weaker interactions and this leads to difficul- 
ties when cooling fermionic vapors. They can be evaded 
for example by sympathetic cooling |ll| with a different 
atom. However it is also feasible to use a scattering res- 
onance, such as a Feshbach resonance, to dramatically 
enhance p-wave scattering. This has been demonstrated 
with 40 K atoms |l2j. The scattering even reaches values 
comparable to s-wave scattering. We will see that this 
means that FQHE gapped states will have characteristic 
energies in the same range as for similar bosonic states. 
At small wavevector, i.e. in the low-energy limit, the 
p-wave phase shift of the two-body scattering problem 
behaves as : 
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where a\ defines the p-wave scattering length. As a con- 
sequence the scattering amplitude is no longer isotropic : 



/i(0) ~ a\k 2 cos6, 



(3) 



where e is the angle between ingoing and outgoing 
wavevectors. We now use an effective potential which 
mimics the behavior in Eq.@ when treated in the Born 
approximation. It is given by : 
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where the quantities r = ri — r 2 and p = |(pi — p 2 ) per- 
tain to the relative particle. The use of such a potential 
is enough for our purpose since we will study the dynam- 
ics of the LLL only : the massive degeneracy is lifted at 
first-order in the potential (higher orders involve Landau 
level mixing). 

We now turn to the quantum Hall regime for fermions 
when u> = u>q in Eq. ■ Assuming a 2D regime with i z = 
^Jfi/mujz the confinement length along z, the interaction 
potential can be written as gf £ 4 p 8^ 2 \r) p where the 
vectors r and p are now 2D with : 
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The coupling constant gf sets the scale of the FQHE 
phenomenon. 

The interaction Hamiltonian in the LLL can be written 
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where m is the relative angular momentum (RAM), 
hence the sum runs over all positive odd integers for 
fermions, and P m (hj) is the projection operator for par- 
ticles i and j onto RAM m. The coefficients V m fully 
characterize the interaction problem in the LLL. They 
are called pseudopotentials after Haldane 0, Q| . With 
the pure p-wave interaction Eq. , only m — 1 scattering 
is allowed and the interaction Hamiltonian reduces im- 
mediately to the pure hard-core model for which V± ^ 
and all the other pseudopotentials are zero. The hard- 
core model ^[ El is known to constitute an 
excellent blueprint of the FQHE. Notably the celebrated 
Laughlin wavefunction |l7| : 
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is an exact zero-energy ground state of the hard-core 
model. In fact this is the most spatially compact zero- 
energy state for the hard-core model. It has a well- 
defined angular momentum L z = 3N(N — l)/2. It cor- 
responds to filling 1/3 of the LLL and thus is the exact 



ground state for ultracold fermions. This is analogous to 
the Bose problem where s-wave scattering leads to the 
Laughlin ground state (with exponent 2) at filling 1/2. 
It is also known that the hard-core model exhibits the 
prominent Jain sequence of incompressible fluids for 
fillings of the form v = p/(2p± 1). We have estimated 
the gaps at some of these fractions by performin g ex act 
diagonalizations in the spherical geometry [ijj, |2(J, Wi . A 
sphere of radius R is threaded by a flux 4irR 2 B which is 
an integral multiple (2S) of the flux quantum by Dirac 
quantization condition. Incompressible fluids appear for 
special matching of the number of particle vs. 2S. We 
obtain the low-lying energy levels for a small number of 
particles and then perform finite-size scaling to obtain 
estimations of the thermodynamic limit . 
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Figure 1: Energy gaps for neutral and charged excitations at 
v — 1/3 , v — 2/5, and neutral gaps for v = 3/7. The lines 
are our best fits and energies are in units of gf. 

For the most stable fluid at v = 1/3 the low- lying neu- 
tral excited states are dominated by a well-defined collec- 
tive density mode and we have obtained the correspond- 
ing gap by studying up to N=13 fermions. The gap can 
also be evaluated from the charged excitations, i.e. the 
Laughlin quasiparticles. For the p-wave problem, there is 
the peculiarity that the quasiholes are gapless at v = 1/3. 
It is thus enough to study the quasielectrons that can be 
nucleated by removing one flux quantum from the ref- 
erence v = 1/3 situation. We find that both estimates 
scale nicely to a common value ~ 0.8 gf : see fig.(l). For 
the fraction v = 2/5 the gap is ~ 0.5 gf obtained by the 
study of systems for N=4...12 fermions. In this case there 
are the two types of charged excitations, quasiholes and 
quasielectrons, each having a nonzero gap and we find 
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that the sum of these gaps converges to a value compati- 
ble with the neutral gap. There are fewer available values 
of the number of particles as we go down the hierarchy 
and for the next fraction c = 3/7 we estimate the gap 
to be ~ 0.4 and for 4/9 it is more difficult to give a 
reliable estimate, the gap is smaller and of the order of 
~ 0.3 gf. For values of the parameters £ z , £ and a\ typi- 
cal of present experiments 0, E3 the § a P s ma Y be of the 
order of the nanoKelvin. 

For fillings less than 1/3, the fermion system does no 
longer display incompressibility because there are prolif- 
erating zero-energy states when we increase the angular 
momentum. Some of these are edge excitations of the 
droplet and for larger angular momentum they are quasi- 
holes. Due to particle-hole symmetry, these very same 
modes fill any gap in the region 1 > v > 2/3. The Jain 
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Figure 2: Energy spectrum for N=9 fermions in the spherical 
geometry at filling v = 1/2. The CFs feel zero flux and can 
be interpreted as forming a closed shell. The inset shows 
the scaling of the gap of the closed shell states N=4,9,16. 
Energies are in units of g and the horizontal axis is total 
angular momentum. 

sequence has an appealing interpretation in terms of com- 
posite fermions (CF). These entities are naively a fermion 
bound to an even number of flux quanta of a fictitious 
field. The total field acting upon the CFs is then the sum 
of the external field and the fictitious field : when treated 
in a mean-field manner this explains the FQHE of elec- 
trons as the integer quantum Hall effect of CFs. The Jain 
sequence has an accumulation point at v = 1/2 : at this 
filling the CF experience zero net flux and form a Fermi 
liquid like ground state |23.l23|. These CFs have remain- 
ing interactions and also an effective mass m* which is 
entirely due to interactions. To estimate this mass, we 
use a special matching of the flux 25 = 2(N — 1) giving 
zero net field on the sphere 0> HH f° r the CFs. In the 
Coulomb case, it has been shown that many features of 



the spectrum can be successfully interpreted by reason- 
ing with free CFs, eventually supplemented by second 
Hund's rule. There are closed shell configurations when 
the number of fermions is a square : N=(£+l) 2 and £ is 
thus the total angular momentum of the highest occu- 
pied orbital. It is the equivalent of the Fermi momentum 
on the sphere and thus we call it £p. We expect that 
the closed shell sequence N=(^f+1) 2 =4, 9,16,. .. display 
ground states with zero total angular momentum and 
should have good scaling properties towards the ther- 
modynamic limit, as is the case for Coulomb interac- 
tions [2^, [2(| . From a closed-shell configuration one can 
form particle-hole (ph) excitations : the lowest-lying such 
excitations is obtained by promoting a fermion from the 
shell with momentum £p to the empty shell at momen- 
tum £p+l, leading to a branch extending from L=l up to 
2£f+1. Above this branch we should find two-particle- 
two-hole states extending up to 4£p and so on. This 
is exactly what we find for the hard-core model. The 
low-lying levels of N=9 fermions are displayed in fig. 2. 
Above the singlet ground state we clearly identify the two 
branches predicted by the free CF model. 

We can obtain an estimate of the CF mass at v = 1/2 
by using the free CF model jl^, |2|| . A free CF on the 
sphere has an energy given by E = 1(1 + l)/(2m*R 2 ) 
where I is the angular momentum. As a consequence, 
the gap of the one ph branch is A = (£p + l)/(m*R 2 ). 
If we fix the density p the scaling law of the gap be- 
comes A N = Air pVN J \m* (N — 1)). In the spherical 
geometry we have to take into account the fact that 
there is a nontrivial shift between the flux 25 and the 
thermodynamic limit value N/v. For finite numbers of 
particles the density is not exactly equal to the thermo- 
dynamic limit value. Better scaling properties |2l| are 
obtained by rescaling the magnetic length by a factor 
\jN/(2Sv) (going to unity for N — > 00). We thus find 
Ajv = Airp/{m*y/N). This scaling is obeyed for the sizes 
N=4,9,16 (see inset of fig. (2)) and this leads to an esti- 
mate of the effective mass : 
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When the number of fermions lie between closed shell 
values we have checked that the ground state angular 
momentum is given by second Hund's rule (maximum 
L), as is the the case for Coulomb interactions [24| . 

Finally we discuss the width of the Hall plateaus in the 
case of trapped atomic vapors. Contrary to the 2D elec- 
tron systems in semiconductor devices there is no source 
of disorder to pin the quasiparticles that are nucleated 
when we deviate from the fine-tuning of a quantum Hall 
fraction. The role of disorder is thus played by the finite 
number of particles of the system and the quantum Hall 
plateaus are expected to be of vanishingly small width in 
the thermodynamic limit. We can give precise estimates 
by considering the equation of state of the rotating sys- 
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tern, i.e. the value of the ground state angular momen- 
tum as a function of the rotation frequency (L z )(lo). At 
the critical frequency, the Hamiltonian is rewritten as a 
magnetic field problem (see Eq.JTJ). If the frequency is 
slightly less, loq — 5u>, then we have the small field —5loL z 
acting upon the purely magnetic problem. It leads to a 
trivial shift of the energies of the FQHE problem that 
will change the ground state when increased. This is 
seen from the typical spectrum displayed in fig. (3) in 
the condition of the critical rotation. The Laughlin state 
at lJf u9hUn = 3N(N - l)/2 is the ground state and the 
first excited state is the quasielectron at lj^ au a hlm — ]\[ 
with a nonzero gap A qe . When adding a —SloL z shift 
the Laughlin will remain the ground state till the quasi- 
electron energy becomes lower for a critical value equal 
to : 
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If we increase 5u beyond this value, quasielectrons are 
nucleated forming a fluid that will condense into a new 
FQHE fluid. We expect the result above in Eq.@ to be 
generic. This picture is essentially dual to the nucleation 
of vortices at small rotation frequency 0. 




Figure 3: Energy spectrum for N=5 fermions in the disk ge- 
ometry as a function of the angular momentum. The Laughlin 
state is the unique zero-energy state at L z — 30. The lowest- 
energy state at L z — 25 is the quasielectron with a finite gap. 

We have shown the appearance of the Jain principal 
sequence of quantum Hall fractions in ultracold rotat- 
ing fermionic vapors. The composite fermion picture 
gives a successful account of the observed fractions as 
well as their collective mode excitations. The gaps we 
estimate from exact diagonalizations are of the order of 
h 2 ai/m£ z £ 4 . At half-filling of the lowest-Landau level, 
there is a Fermi liquid-like state of composite fermions 
and their effective mass m* /m is controlled by £ 2 £ z / (af ). 

We thank Yvan Castin and Jean Dalibard for numer- 



[1 

[2 
[3. 

[4. 
[5 

ir 

[8 
[9 

[10 

[11 

[12 

[13 
[14 

[15 
[16 

li- 
tis 

[19 
[20 
[21 



[22 
[23 

[24 
[25 



Electronic address: Nicolas.Regnault@lpa.ens.fr 
Electronic address: Thierry.Jolicoeur@lpa.ens.fr 
M. R. Matthews , B. P. Anderson, P. C. Haljan , D. S. 
Hall, C. E. Wieman , and E. A. Cornell, Phys. Rev. Lett. 
83, 2498 (1999). 

P. C. Haljan , I. Coddington, P. Engels, and E. A. Cor- 
nell, Phys. Rev. Lett. 87, 210403 (2001). 
K. W. Madison , F. Chevy, W. Wohlleben, and J. Dal- 
ibard, Phys. Rev. Lett. 84, 806 (2000); F. Chevy, K. 
Madison, and J. Dalibard, Phys. Rev. Lett. 85, 2223 
(2000) 

J. R. Abo-Shaeer , C. Raman, J. M. Vogels , and W. 
Ketterle, Science 292, 476 (2001). 

N. K. Wilkin, J. M. F. Gunn, and R. A. Smith, Phys. 
Rev. Lett. 80, 2265 (1998). 

N. R. Cooper and N. K. Wilkin, Phys. Rev. B60, R16279 
(1999). 

N. K. Wilkin and J. M. F. Gunn , Phys. Rev. Lett. 84, 
6 (2000). 

N. R. Cooper, N. K. Wilkin, and J. M. F. Gunn, Phys. 
Rev. Lett. 87, 120405 (2001). 

P. Rosenbusch, D. S. Petrov, S. Sinha, F. Chevy, V. 
Bretin, Y. Castin, G. Shlyapnikov, and J. Dalibard, Phys. 
Rev. Lett. 88, 250403 (2002). 

B. DeMarco, J. L. Bohn, J. P. Burke, Jr., M. Holland, 
and D. S. Jin, Phys. Rev. Lett. 82, 4208 (1999). 

C. J. Myatt, E. A. Burt, R. W. Christ, E. A. Cornell, 
and C. E. Wieman, Phys. Rev. Lett. 78, 586 (1997). 
C. A. Regal, C. Ticknor, J. L. Bohn, and D. S. Jin, Phys. 
Rev. Lett. 90, 053201 (2003). 

F. D. M. Haldane, Phys. Rev. Lett. 51, 605 (1983). 
F. D. M. Haldane, in "The Quantum Hall Effect" edited 
by Prange and Girvin, (1983). 

S. A. Trugman and S. Kivelson, Phys. Rev. B31, 5280 
(1985). 

A. H. MacDonald, in "Mesoscopic quantum physics", Les 
Houches, E. Akkermans, G. Montambaux, J. L. Pichard 
editors, Elsevier Science, Amsterdam (1994). 
R. B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983). 
J. K. Jain, Phys. Rev. Lett. 63, 199 (1989). 

F. D. M. Haldane and E. H. Rezayi, Phys. Rev. Lett. 54, 
237 (1985). 

G. Fano, F. Ortolani, and E. Colombo, Phys. Rev. B34, 
2670 (1986). 

R. H. Morf, N. d'Ambrumenil, and S. Das Sarma, Phys. 
Rev. B66, 075408 (2002); N. d'Ambrumenil and R. Morf, 
Phys. Rev. B40, 6108 (1989). 

A. Lopez and E. Fradkin, Phys. Rev. B44, 5246 (1991). 

B. I. Halperin, P. A. Lee and N. Read, Phys. Rev. B47, 
7312 (1993). 

E. Rezayi and N. Read, Phys. Rev. Lett. 72, 900 (1994). 
R. Morf and N. d'Ambrumenil, Phys. Rev. Lett. 74, 5116 
(1995). 

[26] M. Onoda, T. Mizusaki, T. Otsuka, and H. Aoki, Phys. 
Rev. Lett. 84, 3942 (2000). 



